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Abstract 

The XXZ Heisenberg chain is considered for two specific limits of the anisotropy pa- 
rameter: A —7- and A — > — oo. The corresponding wave functions are expressed by 
means of the symmetric Schur functions. Certain expectation values and thermal cor- 
relation functions of the ferromagnetic string operators are calculated over the base of 
A^-particle Bethe states. The thermal correlator of the ferromagnetic string is expressed 
through the generating function of the lattice paths of random walks of vicious walkers. 
A relationship between the expectation values obtained and the generating functions 
of strict plane partitions in a box is discussed. Asymptotic estimate of the thermal 
correlator of the ferromagnetic string is obtained in the limit of zero temperature. It 
is shown that its amplitude is related to the number of plane partitions. 
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1 Introduction 



1.1 XXZ Heisenberg chain 

A system of spin 1/2 particles occupying sites of one- dimensional lattice, widely known as the 
quantum XYZ Heisenberg chain [TJ, has attracted considerable attention both in theoretical 
and mathematical physics, and it has been thoroughly investigated for a long time [2H8]. 
The Quantum Inverse Scattering Method, developed for solving the integrable models of 
quantum field theory and statistical physics [HI EE], has been also used to investigate the 
XYZ Heisenberg chain [ll,12j. An important special case of the XYZ model, so-called, 
XXZ spin chain also attracts considerable attention [T3 | fl3 |fT8l - [20] . The Hamiltonian of 
the XXZ magnet has been diagonalized by the coordinate Bethe Ansatz Method in [3H5]. 
The Algebraic Bethe Ansatz has been used in [11] to solve the XXZ model. The problem 
of calculation of the correlation functions of the model in question in the framework of the 
Algebraic Bethe Ansatz has required serious efforts [T5HT7] . 

The random walks is a classical problem both for combinatorics and statistical physics. 
The problem of enumeration of the paths made by the, so-called, vicious walkers on the 
one- dimensional lattice has been formulated and investigated in details by Fisher [21]. The 
problem mentioned still continues to attract considerable attention both of physicists and 
mathematicians [22H33]. 

Random walks on one-dimensional periodic lattice can be related to the correlation func- 
tions of the XX Heisenberg magnet. Certain operator averages taken over the ferromagnetic 
state of the XX model play a role of the generating functions of the number of paths traced 
by the vicious walkers. The problem of enumeration of trajectories of vicious walkers by 
means of the correlation functions of the XX model has been studied in the series of pa- 
pers J35H3B]- The approach of Refs. [5oT[5%] will be further explored in the present paper. 
We shall consider two limits of the XXZ spin chain from the point of view of random walks 
of vicious walkers, as well as from a viewpoint of enumeration of boxed plane partitions. 

Let us begin with the XXZ model defined on one- dimensional lattice consisting of M+l 
sites labeled by elements of the set M = {0 < k < M, k e Z}, M + 1 = (mod 2). The 
corresponding spin Hamiltonian is defined, in absence of external magnetic field, as follows: 

1 M A 
Hxxz = -- + + y K+i° Z k ~ 1)) , (1) 

k=0 

where the parameter A 6 1 describes the internal anisotropy of the model. For instance, 
the choice A = ±1 corresponds to, so-called, isotropic XXX spin chain solved in [2J. The 
local spin operators a h = ~(af ± ia v k ) and c|, dependent on the lattice argument k E Ai, 
are defined as (M + l)-fold tensor products as follows: 

Crf = CT°<g>---®CT°®^cr^®cr <g>---®cr , (2) 

k 

where cr° is 2 x 2 unit matrix, and <r* at k th place denotes a Pauli matrix, a* G su(2) 
(superscript # implies either x,y, z or ±). Therefore, the spin operators act over the state- 
space $)m+i given by the tensor product of M + 1 copies of the linear spaces l)k = C 2 : 

M 

$)ai+i — (8) The commutation rules for the spin operators are given by the relations: 

fc=0 

[ <rt, <h \ = 5 M °i > [ < } = ±2 5 k ,i of . 
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The linear space C 2 is spanned over the spin "up" and "down" states (|f) and respec- 
tively) providing a natural basis so that 

1 \ 1 1 \ / 



M 

The space $)m+i is spanned over of the state- vectors ® |s)fc , where s implies either | or |. 

fc=0 

The periodic boundary conditions crf + r M+1 \ = erf are imposed. 

To represent iV-particle state- vectors of the model, \ ^n(u\, . . . ,Un)), let the sites with 
spin "down" states be labeled by the coordinates ^ 1 < i < N. These coordinates form 
a strict partition /j, = (yU 1; /i 2 , • • • ,^n), where M > ^ > fi 2 > ... > /ijv > 0. There is a 
correspondence between each partition and an appropriate sequence of zeros and unities of 
the form: {e^ = e k(^)} keM , where = Sk )fln , 1 < n < N. This correspondence enables 
another convenient expression for the strict partitions: /x = (M eM , . . . , l ei , e °). It is meant 
here that any non-negative integer k e M. appears times in the present configuration, and 

M 

the condition e k = N is respected. The Hamiltonian ([1]) is diagonalized via the ansatz: 

k=0 

M 

|vmu)>= Yl x^nu^rm, (3) 

where the sum is taken over C^- strict partitions /x. The state Iff - ) in ([3]) is the fully polarized 

M 

state with all spins "up": |f|~) = (^) |t)n- Besides, it is proposed to use bold-faced letters 

n=0 

as short-hand notations for appropriate iV-tuples of numbers: for instance, u instead of 
(tii, • • • ,%), etc. Therefore, the wave function x XXZ ( u ) i n © is of the form: 

x xxz n = £ ^(u)<^f...^7, (4) 

where summation goes over all elements of the symmetric group of permutations S PIiP2 PN = 
1. 2 



s( ' ' '"' V The amplitude ^.5 is given by the product: 
Vpi, r»2, • • • i PnJ 



^ u = n — 2 2 pa pj ■ (5) 

X X nt* — qt* 

l<j<i<N Pi Pi 

The state- vectors (jSJ) are the eigen-states of the Hamiltonian (JTJ, 

ffxxz |^iv(u)) = ^(u) |^jv(u)), 
if and only if the variables u\ (1 < I < N) satisfy the Bethe equations 

The corresponding eigen-energies are given by 

En(u) = --$> 2 + M 7 2 -2A). (7) 
i=i 
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1.2 Outline of the limiting models and the problem 

In our paper we shall consider two special cases of the XX Z model, namely the A — > and 
A — > — oo limits. The A — > limit known as the XX Heisenberg chain is a most popular 
and studied one [39ti4Tj . The Hamiltonian of the XX model describes the nearest-neighbor 
interactions of spin "up" |t) and spin "down" \l) states located on sites of the periodic chain 
in zero magnetic field as follows: 

1 M 

H XX = - ^ LMk+i°k + a t+i a k) ■ ( 8 ) 

fc=0 

It is crucial that the system described by the Hamiltonian ([8]) is equivalent to free fermions 
[4"T] . The case of the XX magnet can also be deduced by taking the limit of infinite on- 
site repulsion in the boson Hubbard model (the "hard-core" bosons) [48] . Therefore, the 
XX model is interesting for investigating the (quantum) phase diagram of the Hubbard 
model, as well as for description of the Frenkel excitons jl9]. In the last few years, the XX 
model has attracted attention in connection with the quantum information and computation 
theory [SHEE]. 

The state-vector of the XX model and the correspondent Bethe equations are obtained 
at A = from (J3J) and respectively Up to an irrelevant pre-factor, the wave functions 
of the model are equal to 



and the Bethe equations are 



Kn<KN 



uf M+1) = (-l) N -\ l<j<N. (10) 
The substitution u 2 = e j brings these equations to the exponential form, 

e m + i )ej = (_!)"-!, (ii) 

with the solutions: 

J,= (12) 

3 M + 1V 3 2 /' y 1 

where Ij are integers or half-integers depending on whether iV is odd or even. It is sufficient 
to consider a set of iV different numbers Ij satisfying the condition: M > Ii > I 2 > ■ ■ ■ > 
In > 0. The notation 6 for iV-tuple (9\, 9 2 , . . . , 6 N ) of solutions (TT2"j) will be especially 
convenient for usage below in order to stress that one is concerned with the solution of the 
Bethe equation. Otherwise, it is appropriate to use u as an indication that arbitrary set of 
parameters is meant. It follows from (0) that the eigen-energies of the XX model are equal 
to 

E**W = -f>9, = -£00.(^(1,-^)). (13) 
i=i i=i v J 

The ground state of the model corresponds to the following solutions of the Bethe equations: 

0i = —(N-3-—), l< j < N. (14) 
3 M+ 1 V J 2 J ~ K ' 
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Less studied limit of the XXZ model is the Strong Anisotropy (SA) limit A — > — oo 
[UEM22H52] ■ I n this limit the behavior of the system is described by the effective Hamiltonian 
which is formally equivalent to the XX Hamiltonian supplied with requirement forbidding 
two spin "down" states to occupy any pair of nearest-neighboring sites [52|l53] : 



fc=0 

where the projector V cuts out the states with the spin "down" states at any pair of nearest- 

M 

neighboring sites: V = J^[(l — Qk+iQk)- The local projectors onto the spin "up" and "down" 

fc=0 

states are equal to: 

Qk = \ (4 + 0-1), % = \ (4 -a z k ), q k + q k = I, k G Ai , (16) 

where the operators af are defined by (J2J). 

In the limit A — > — oo, the wave function (J5J) takes the form: 

Xl\u) = det(uf^- N+ %^ k < N J] ^i-<T\ (17) 

l<n<K7V 

where the coordinates of the spin "down" states form strict decreasing partition fx (as in (jj]) 
and (jHD), i-e., M > \i\ > /x 2 > • • • > Hn > 0. It follows from ( TT7|) that the wave function is 
not equal to zero if and only if the elements 1 < i < N, satisfy the exclusion condition: 
fii > + 1. It is crucial that in the considered limit the occupation of nearest sites is 
forbidden, and the hard-core diameter equal to duplicated inter-site separation arises. The 
Bethe equations of the model take the form: 

N 

e i{M+i-N)e k = Y[ e - i0 i , 1 < k < N , (18) 

and have the solutions 

3 M+l-N\ 3 2 /' K J 

1 N 

where P = — 9j, and Ij are integers or half-integers depending on N being odd or even, 

n 3=1 

satisfying the condition M — A^>/i>/ 2 >--->^w>0. The ground state of the model is 
defined by the solutions 

while P = is fulfilled. The eigen-energy of the model is: 

N / \ 

^-Ecos(-i^^-^-P)). (2D 
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Similarity between these two limits, (JHJ and f lT5|) . is that their wave functions are ex- 
pressible through the Schur functions [56J : 



<? M ~ 9 (r t t ) ~ ( j ll^^L 
= o\{Xi, x 2 , ■ ■ ■ , x N ) = — 7777~jv-fc 



det(Xj )i<j,k<N ^ 
= det(x^ +Ar ~ fc )i< J - fc <iv f { 

l<n<J<iV 

where A denotes the partition (Ai, A 2 , . . . , A at) being iV-tuple of non- increasing non-negative 
integers: £ > Ai > A 2 > ... > Ajv > 0. Indeed, any strict partition M > \i\ > fi 2 > ... > 
> and non-strict partition M + l — A^>Ai>A2>--->Aat>0 (denoted as fx and 
A, respectively) can be related by means of the relation Xj = fij — N + j , where 1 < j < N. 
In other terms, A = fi — d, where 8 is the strict partition (N — 1, iV — 2, . . . , 1, 0). So the 
wave function of the XX model (JSJ may be represented as 

X^(u) = S A (u 2 ). (23) 

Any strict partition /x with the elements respecting the exclusion condition /ij > yU i+1 + 1 is 
connected with the non-strict partition A by the relation A = fi — 2S, where M + 2(l — N) > 
Ai > A 2 > ■ • ■ > Aat > 0. Therefore, Eq. ( Tl7|) is re-expressed: 



= %(u 2 ) . (24) 

It is useful to remind a graphical picture (see Figure 1) of the correspondence between 
strict partitions /x and non-strict partitions A used in (12"3"|) . Namely, to each partition A we 
associate a set of iV-tuples (3(A) as follows 



0(A) = {\ j -j + l\i<j<N}cZ+ \ 



2 

On another hand, any non-strict partition A can be represented as a rectangular table (the 
Young table) consisting of iV columns so that Aj, Vi, is the height of i th column (Aj < 
M-N + l). We shift each element of the set 0(A) by JV + \. Then we assign the numbers 
thus obtained to the projections along the vertical dashed lines onto the horizontal axis. 
The set of points on the horizontal axis just provides the strict partition ii. For instance, 
the diagram on Figure 1 is drown for M = 8 and N = 4, and we have got respectively: 
A = (5,3,2,2) and [i = (8,5,3,2). 

In our paper we shall study the thermal correlation function of the states (would be called 
as ferromagnetic strings) with no spins down on the last n + 1 sites of the lattice. We are 
going to consider the expectation value defined by the ratio: 

TIB R\ - {9N(e)\JKe-^JK hMg)) ft _ TT - ton 

™ = — <&m \*m — • n " = 3J' () 

where G C, and the projector Tl n is expressed by means of qj ffTBT) . Besides, H in T(6, n, (5) 
fT25|) implies either Hxx © or Hsa (|T5l) . and 6 indicates that the eigen-state | ^n(0)) is 
calculated for solution of the Bethe equation fTTTj) or (TTBl . respectively. Our calculations will 
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Figure 1: The strict partition // 2 , fa, fa) and the corresponding Young table. 



be based on the similarity of the problem under consideration with the problem of enumer- 
ation of admissible lattice paths made by random vicious walkers. We shall extensively use 
the relation for the Schur functions which is due to the Binet-Cauchy formula [57] : 

S x (x 2 1 ,...,x 2 N )S x (y 2 1 ,...,y 2 N ) (26) 

XC{L N } 

= det(T jfc )i<j,fc<iv n {vj-viy 1 n k - 1 • 

±<k<j<N l<m<l<N 

The entries of the matrix Tjk take the form: 

1 - (x m r^ 

• (27) 

Summation in (126 p goes over all non-strict partitions A into at most N parts so that each is 
less than L: L > Ai > A2 > ■ ■ ■ > Xn > 0. The notation for the range of summation in fl26|) 
will be extensively used in the rest of the paper. 



2 XX Heisenberg chain 

2.1 The Bethe states and form-factors 

Before to proceed with calculation of T(0,n,f3) (|25|) . let us apply our approach to more 
familiar examples. First of all, let us specialize, with regard at (123|) . the state- vector fl3]) and 
its conjugated as follows: 



i^(u))= e s A (u 2 ) n K) efc i^> , 

M ^ 

Nv)i= e (tin^^v- 2 ), 

AC{(A/+l-Af) JV } k=0 



where summation goes over all non-strict partitions A, which are related to the non-strict 
partitions /j, = A + 8. Because of the orthogonality relation 



M M M 

6 n 6 n 1 



fc=0 1=0 n=0 

the scalar product of the state-vectors (1281) takes the form: 

W1H»H)= £ s^)s^) = d ;* ( g;^: j r . (29) 

AC{(M+l-Af) JV } ' 

We use in (129]) the notation for the Vandermonde determinant, 

V(u 2 ) = J] (u 2 - 1&) , (30) 

l<m<Z<7V 

while the entries T^j of the matrix of the size N x N are of the form: 

Equations ( 121?]) and fl3"TT) are specifications of Eqs. f[2"6"]) and f[2"T]) . respectively. On the 
solutions ( TT2|) . the entries ( |3T|) are equal to Tj^ = (M + l) A jfc, where l'Hospital rule is taken 
into account on the principle diagonal. Let us use the exponential parametrization for the 
solutions of the Bethe equations in compact form: 

u 2 = e w , e ie = (e idl , e^ 2 , . . . , e i0N ) , (32) 

where the "angle" notation ( Fl2|) or (fl9|) is meant. Thus, the answer for squared norm 
A/" 2 (0) ee (ty N (0) \ ^n(0)) of the Bethe eigen-vector ([28]) arises as follows: 

Af*(0) = {M + 1)N = {M + 1)N (K) 

1 ] V(e»)V(e-*) EI 2(1 - cos(0i - m )) ' 1 J 

l<m<«<7V 

Let us now turn to the ratio: 

T(V ' U '" ) = AT(v)AA(u) ' (34) 

where the projector Tl n is defined in and 7V 2 (u) = (^/^(u) | \1/jv(u)) for arbitrary 

parametrization of the state- vectors. With regard at fT2"8l) . we calculate: 

M M-n-l 

U n |<Mu)> = £ 5 A (u 2 ) ( J] (<r")°) ( n K)^) lit), (35) 

AC{(Af-Af-n) JV } k=M-n k=0 

where summation goes over non-strict partitions A respecting the condition: M — N — n > 
Ai > A2 > ... > Aat > 0. Taking into account (128]) and (|35|) . and using the Binet-Cauchy 
formula, we calculate the nominator of ( 134]) : 

(vMv) I 77 n I <Mu)) = ^ 5 A (v- 2 )5 A (u 2 ) 

AC{(M-AT-n) JV } 



1 , ,7 1 - (&M) M _z 

l<j,k<N 



det ^^7^ ■ (36) 



V(u 2 )V(v- 2 ) ^ 1 - u\jv) 
8 



Assume that the sets of the parameters v and u in (134]) coincide and consist of the so- 
lutions (TT2"]) . Then, expression T(0,n) = T(e 40 / 2 , e* 0//2 , n) is related to, so-called, Emptiness 
Formation Probability, which provides the probability of formation of a string of the empty 
(i.e., spin "up") states on last n + 1 sites of the lattice. Eventually, Eq. (136]) leads to T{6, n) 
in the same determinantal form as in [1T1H3] (see therein for more Refs.): 

ne.n) = d*((l - — )fc + (M+1)(1 _ e , ( ^ j)) d - W)^, (37) 
where the parameters 6i, 1 < I < N, correspond to the parametrization fTT2|) . 

2.2 Thermal correlator of ferromagnetic string and random walks 
of vicious walkers 

Let us turn to calculation of the following expectation value: 

T(v ' u ' n ' A = A7(vWH ' (38) 

which is clearly reduced to 7~(v, u, n) (134)) at /3 = 0. However, our aim is to obtain T(0, n, (3) 
( 125]) just taking the arguments u and v in (138]) coinciding with the same solution (fl2|) . 

Let us use the technique presented above to calculation of the nominator of (138]) . Using 
(13"5"]) (and its conjugated), one gets: 

(*jv(v) | n n e-^ xx n n \ V N {u)} 

(39) 

E S A ,(v- 2 )S Afl (u 2 )F^(/3). 

A L ,A fl C{(Af-V-n) JV } 

Summations in (139|) run over non-strict partitions \ L and X R of the same kind as in (1351) . Su- 
perscripts L and i? are only to distinguish two independent summations. The corresponding 
strict partitions ii L and fj, R are defined as follows: /j, l ' r = X L ' R +6, where d = (Si, 6~ 2 , ■ ■ ■ , 5jv), 
Sj = N — j . The notation F^l.^r^) implies the following average: 

(40) 

which is nothing but 2iV-point correlation function over the ferromagnetic state. It is related 
to enumeration of admissible trajectories which are traced by N vicious walkers traveling 
over sites of one- dimensional chain [351438] . 

Indeed, let |P^(/if , . . . , n R — > . . . , //^)| be a number of trajectories consisting of 
K links made by iV vicious walkers in the random turns model. Here, the initial and final 
positions of the walkers on the sites are given respectively by elements of the strict decreasing 
partitions /if > \i R > ■ ■ ■ > fi R and \i\ > ^ > * ' ' > ^n- We introduce the notation Vf for 
the operator of differentiation of K th order with respect to i at the point i = [35] • Then, 
the application of the operator T>^ 2 to the correlator (T4T3]) results in the average of the type 

<%VS • • • (-2^*)^. V • • • o-^lt) ■ 
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This average provides the numbers \Pk{^ r , • • • , A*jv — lAi ■ ■ ■ ■• ^n) I > as ^ has been established 
in [36] with the help of the commutation relation 



N 



iHxx, . . . a^] = £ %f . . . [H xx , ^\a^ . . . . (41) 



fc=i 



The condition of non-intersection of trajectories of the walkers is expressed by the vanishing 
of the correlation function ( 140]) for any pair of coinciding indices y R or . Thus we conclude 
that the average (1591) turns out to be the generating function of the polynomials dependent 
on 2N variables, u\,u^ ■ ■ ■ , u 2 N and v f 2 , v^ 2 , • • • , v^ 2 , as follows: 



(*jv(v) \n n e-f* 6 xxn n \* N {u))^ 

(42) 

E |Px(^^^)|5 A ,(v- 2 )^(u 2 ) 

A L , \ R C{(M-N-n) N } 

(remind that fi L,R = X L ' R + S). As it is shown in [3B], the number of trajectories consisting of 
K links, which are traced by iV vicious walkers on an axis, i.e, \Pk(h — > A* L )|, is expressed 
through the number of trajectories of the same "length" K, which are traced by a single 
walker traveling over sites of iV-dimensional lattice of infinite extension. 
The correlator ( I40p respects the following equation: 

^m^,l4,...,fi%;fj,f,^,... >f t§(^) = 

l n (43) 

= — > (F..L..L ,,L.,,R,,R „Sii ,,r(B) + F,,L ,,L ..L...R..R ,,R i ,,r(B)). 

fe=l 

Equation ( 143|) has been considered in [36] for the case of the periodic boundary condition 
with respect to the lattice argument and with the initial condition: 

N 



k=l 



Solution to (l4"3"j) can be expressed as the determinant of the matrix { F k;i{P)) 1<kl<N [35|l36]: 

F ti,tf,...,fotf,tf,...,,*%(P) = det i F rt;tf(P))i< k ,i<N> ( 44 ) 
where the entries respect the following difference-differential equation: 

^ F M = \{P k+ iAP) + F k-iM)- (45) 

Similar equation can be also obtained for the fixed index I. 

It can be checked that the transition amplitude (ff \a k e~@ Hxx G^\ f|~) respects (l4"5j) [351438] . 
This average can be considered as the generating function of the number of walks with 
random turns of a single pedestrian traveling between I th and k th sites of (periodic) chain 
[36|,[38]. Solution to (j4"5l) can be written as the following sum: 

M 

F k-,i{0) = -—^eP™** e^ k - l \ (46) 
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2tt / M\ 

where the parametrization (f) s = — — - ys — —j is used. The periodicity condition with 

respect of the lattice argument and the "initial" condition Fj..i(Q) = 5k j are imposed. 

Using (jlijl and fj46l) . we re-express (JlOj) through the Schur functions fl22|) and the Van- 
dermonde determinants fl30|) as follows [35] : 

JV 

^ 1 ^-^ /3E cos (4,) 

Fu?"u.*{0) = 777 777 > e i=1 

,M v 1 (M+1) N ^ (A^\ 

xV(e^)V(e-^) S' x z.(e i *)5 , A «(e- < *) . 

The parametrization S , < s < M, used in fH7|) is the same as in ( T4T)j) . We continue 
to use bold-faced letters to denote iV-tuples of numbers: for instance, </> corresponds to 
(0fen 0fc 2 5 ■ ■ ■ 5 0fciv)- We substitute ( H7|) into (13"9"1) . and use flHET) in order to calculate the sums: 

P(v- 2 ,e^) = ^5 Ai (v- 2 )5 A ,(e^), V{e^, u 2 ) = ^ 5 A «(e^) 5 A «(u 2 ) . (48) 

A^ A« 

The range of summation in (j4"81) is taken as follows: A L , A R C {(M — iV — n) - ^}. Then we 
obtain: 

(^(v)|i7 n e-^i7 n |^v(u)> = 

N 

1 «, J2 cos(0 fc; ) 

2^ e 1=1 x (49) 



(M + l) 7V V(u 2 )V(v- 2 / 

v / V / V / M>k 1 >k 2 - >k N >0 

1 - (e^^7 2 ) M -"x /l - (^e-HfcjW-n 

x det r- 1 — a — det ' 



1 — e^Vj 2 J i<i,j<N \ 1 — %(t>k i /i<p,i<N 
The Binet-Cauchy formula enables to evaluate ( I49p as follows: 

(^jv(v) |77 n e-^i7 n |^(u)) = 

x fc ''- u 7 l<i,j<N 

where F k .j(f3) is defined by ( 145]) . Clearly, the relation (15"U|) at /3 = is reduced to (1351) . 
Expression for T(0, n, /3) (125 p can straightforwardly be obtained from (1501) at coinciding 
parameters u = v (being solutions of the Bethe equations) with the help of M 2 ( l33~j) : 

1 /M-n-l 

T(0,n,P) = det £ F M (/3)e^"^) 

^ ' \ k,l=0 / l<i,j<N 

3 Strongly anisotropic XXZ chain 

3.1 The Bethe state vectors and their form- factors 

Let us turn to the strong anisotropy limit, A — > — oo, which is described by the Hamiltonian 
Hsa (fI5j) . The corresponding state- vector is given by Eq. 

M 

\*n(u)) = Yl ^(u 2 ) nK) efe It) , (51) 

AC{(A/-2(Af-l)) JV } fc =0 
11 



where S^(u 2 ) is given by fl24|) . and summation goes over all non-strict partitions. Summation 
over strict partitions \x with the elements respecting the condition /ij > + 1 is equivalent 
to that over the non-strict partitions A, where A = fx — 2d and M + 2(1 — N) > Ai > A2 > 
■ ■ ■ > X N > 0. The scalar product of the state- vector (I5T1) to its conjugated (defined similarly 
to ( 12"8|) ) is given by the relation looking, in turn, similarly to 



\C{(M-2(N-1)) N } 

where the entries Tjy take the form: 



1 - (ul/v]) 1 

J-kj — 1 2/2 ' (.""J 

and the notation fl30l) for the Vandermonde determinant is used. When the sets of parameters 
v and u in left-hand side of fl52]) coincide, usage of the Bethe equations f TT8|) enables to express 
the entries ( |53|) as follows: Tjfc = 1 + (M — N + l)5jk- Therefore the norm of the Bethe 
eigen- vectors Af 2 (9) = (^ N (9) \ ^ N (0)) is given by 

(M + 1)(M + 1 - N) N -' {M + l){M+l-N) N ^ 
lJ V(e*)V(e-*) EI 2(l-cos(^-^ m ))- 1 J 

l<m<i<iV 

The exponential parametrization ( Fl9|) is meant in ( 15^|) in the compact form ( |32|) . It can 
be shown that the scalar product (^^(v) | ^^(u)) (|52|) vanishes (i.e., the state-vectors are 
orthogonal) provided the parameters u and v are independent Bethe solutions. 

The nominator of the ratio (I34p is calculated in the same way as in the Section 2.1, i.e., 
by means of the Binet-Cauchy formula: 

OMv) I n n 1 <Mu)> = J2 %(v- 2 )%(u 2 ) 

XC{(M-2N-n+l) N } 

" V(U 2 )V(V- 2 ) \ I -Ul/v] Jl<j,k<N- [ } 

After use of (1541 and fl55l) . the answer for fl34l) . taken on the solutions of the Bethe equation 
(TT9]) (i.e., the Emptiness Formation Probability), appears in the following form: 

, M-N + l , /, n N P 

+ ( M -iV + l)(l-e^-^)) (1 " 5jfc) ) i<fej<iv ' (56) 
3.2 Thermal correlator of ferromagnetic string 

Let us turn to obtaining of the average T(v, u, n, f$) fl38|) . where the Hamiltonian Hxx is 
replaced by H$a (fl5i) . Using floTl) we obtain the corresponding nominator of fl38|) : 



iy N (v)\n n e-^n n \y N (u)) = 

E ^(v- 2 )%4u 2 )^ ; ^(/3), 

A L , X R C{(M-2N-n+l) N } 
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where F^l.jir(^) is 2iV-point correlation function over the ferromagnetic state given, prac- 
tically, by ( HOjl excepting that Hxx dH]) is replaced by (1T51) . Summations in (1571) go 
over non-strict partitions X L and \ R of the same kind like in (|52|) . The corresponding strict 
partitions J1 L and J1 R are defined as follows: jl L,R = X L ' R + 2S, where d implies the partition 
(N — 1, iV — 2, . . . , 0). It is crucial that now the lattice indices J1 R respect the exclusion 
requirement: occupation of nearest sites is forbidden. Besides, an analogue of the relation 
E2l) can be written as follows: 

£ \p K (ji R ^fi L )\s~ xL (v^)s~ xR (u 2 ). (5t ' 

\ L ,X R C{(M-2N-n+l) N } 

The solutions of the Bethe equation (ITgj) constitute a complete set of the eigen-states [53] . 
Taking into account the orthogonality of the corresponding states, one can consider the 
resolution of the identity operator: 

I = ^AT- 2 (0) \* N (0))(* N (0)\, (59) 
W 

where summation goes over all independent solutions of the Bethe equation ()18|) . and the 
square of the norm Af 2 (0) is given by (15^1) . We shall calculate ( 15T|) inserting (159]) into 
left-hand side of (15T1) and using appropriately Eq. (155|) . We take into account that 



(^v(v) | e"^ SA \* N {0)) = (^v(v) \V N (0))e 



where the expression for the energy En(6) is given by the relation (12T1) . Further, we use 
(1551) and obtain: 



(^v(v) l^e-^A^i^u)) 



1 



(M + 1)(M + 1 - iV)^- 1 V(u 2 )V(v- 2 ) 



e 



-PE N (0) 



M-N>h>l 2 ->l N >0 (60) 



x det ( 3 —r — * ) det ' 1 1 



l — e idt v, 2 J i<i,j<N V 1 — ufe^ e p A<Z,p<AT 

J L 

where summation goes over the ordered sets {Ik}i<k<N, which parametrize the solution 
(lT9~j) . Expression for T(0,n,f3) ([UJ can be obtained by means of (15^1) and (1501) . where it is 
necessary to put u 2 = v 2 = e 10 , while 6 is the solution (120]) of the Bethe equation for the 
ground state. 



4 Boxed plane partitions 

We shall show that the scalar products of the state vectors and the emptiness formation 
probability are related to the generating functions of the boxed plane partitions. 

An array (^i,j)i,j>i of non- negative integers that are non- increasing as functions of both i 
and j G {1, 2, . . . } is called a plane partition iz [S5J. The integers 71"^ are called the parts 
of the plane partition, and | tt | = -k^ is its volume. Each plane partition has a three 
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dimensional diagram which can be interpreted as a stacks of unit cubes (three-dimensional 
Young diagram). The height of a stack with coordinates is equal to the part of plane 
partition 7Tjj. If we have i < r, j < s and 7Tj j < t for all cubes of the plane partition, it is 
said that the plane partition is contained in a box with side lengths r,s,t. If TTjj > 
i.e. if the parts of plane partition 7r are decaying along each column, then 7r is called column 
strict plane partition. We shall call the partition 7r that are decaying along each column 
and each raw (jr^ > 7r i+1 j and 7Tjj > vrj as the strict plane partition. The element tti .1 
of the strict plane partition 7r satisfies the condition 7Ti 1 > 2r — 2, if all 2, j < r. 

An arbitrary plane partition in a box r x r x t may be transferred into a column strict 
plane partition in a box r x r x (t + r — 1) by adding to an array (fti,j)i,j>i the r x r matrix 



7T 



cspp 



/ r 

r 

V 



1 

2 



1 

2 















1 \ 

2 

y 



which corresponds to a minimal column strict plane partition. The volumes of the column 
strict plane partition and correspondent plane partition are related 



7T 



cspp I 



k + 



N 2 (N-1) 



(61) 



An arbitrary plane partition in a box rxrxt may be transferred into a strict plane partition 
in a box r x r x (t + 2r — 2) by adding to an array >i the r x r matrix 



7T 



spp 



/ 2r-2 2r-3 
2r-3 2r-4 

V r — 1 r — 2 



r - 1 \ 

r-2 







which corresponds to a minimal strict plane partition. The volumes of the strict plane 
partition and correspondent plane partition are related 



7T 



spp I 



|tt| +N 2 (N- 1). 



(62) 



The partition function of the three dimensional Young diagrams, or saying differently the 
generating function of plane partitions is equal to 



{*■} 



(63) 



where q is a weight, and summation is performed over all plane partitions in a box. Formulas 
(16T|) and fl62|) provide the connection between partition functions of the plane partitions of 
different types. The generating functions of the column strict and strict plane partitions 
placed into a box N x N x M are equal respectively to 



y cspp 



^spp(g) 



M 2 {N- 



-» n 



l<7,fc<iV 



1 _ gM+i+j-k 
1 _ q j+k-i 



l<j,k<N 



-i) tt 1 - g A w- fc 
11 1 _ gj+k-i ■ 



(64) 
(65) 
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The scalar product fl29|) is related to the partition function of the column strict three- 
dimensional Young diagrams placed into N x N x M box. Really, the parametrizations 
Vj = q~i and Uj = give: 



(* N (q- 5, . . . , q-T) | y N (l, . . . , q— )) = 

E 5 A (g,...,^)5 A (l,...,g JV - 1 ) = (6g) 



\C{(M+1-N) N } 



1 , /l-sJ +fc - 1 



det 



(- 



V(q, q N )V(l, q N ~ 1 ) VI - q j+k ~ x J i<j,k<N ' 
where s = q M+1 , and 

V- 1 (g,... ) ^)V- 1 (l,...,^- 1 ) = g-^- 1 ^- 1 ) II (l-^-*)" 2 - (67) 

l<k<j<N 

The determinant in fl66|) was calculated in the paper [61 J in connection with the alternating 
sign matrices enumeration problem: 



detf 



1 - s^- 1 



\1 — i<j,k<N 



N 



J 

l<k<j<N k,j=l * 



(6* 



Taking into account (l64j). flBTj) . and f l68|) . we obtain for (job]) : 

(* N (q-*, ...,q~) | ^(1, • • • , g— )) = g-- {JV - 1) Z CS pp(g) . (69) 

Thus, Eq. flOSJ) reads that the scalar product of two state-vectors coincides at q — 1 with 
the number of column strict partitions in a box N x N x M, i.e., with Z cspp (l). 

The same parametrizations Vj = q~z and w,- = q^~ enable to express the scalar product 
f)52p corresponding to the strong anisotropy limit. The same representation is valid though 
s = q M ~ N + 2 now, while the range of summation over A takes the form: 

A C {(M — 2(N — l)) N } . 

Taking into account Eq. (|68p . we now obtain: 

<<Mg" V ■ ■ , ) I ^(1, . . . , g^ 1 )) = ff-^^^ppC?) , 

where Z spp (g) is the generating function of strict plane partitions ( j65|) . The corresponding 
value of the scalar product at g = 1 coincides with the number of strict plane partitions 

^spp(l). 

Now lets turn to the expectation value of the ferromagnetic string (136|) . In the present 
parametrization, we obtain (with regard at (IB"?]) and (16"B"]) ): 



(*Ar(g _ 5, . . . , q -f ) | i7 n | ^(1, . . . , g^r)) = 

E S x ( q ,...,q N )S x (l,...,q N -l) = I [ ; ~ ^ 



AC{(A/-JV-n)"} fcJ=1 1 - Q J 
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where s = q M ~ n , i.e, flTDj) differs from fl66l) in the sense that n + 1 is subtracted from M + 1. 
The box containing the plane partitions is now of smaller maximal height: N x N x (M — 
Ti—l). The expectation value in left-hand side of f!7U|) ). being considered at q — > 1, "counts" 
the number of plane partitions of smaller height. The same is true for the case of strong 
anisotropy, and (155|) takes the form: 



CM? - *, • • • , I I *jv(1, ...,gT-)) 



E 5x(?.---.g Ar )5'x(i,-,^- 1 ) = II T 

AC{(M-27V-n+l)^} fc J= l 1 



* (71) 



sg 



_ gj+k-i 



where s = q 



M—N-n+l 



5 Low temperature asymptotics 

Let us go over to the case of long enough XX-chain, i.e., M ^> 1 while N is moderate: 
M ^> N ^> 1. Now the correlator Fk ; i(f3) is approximately given by the modified Bessel 
function instead of (1461): 



7T 



= = ^- f e^e^-Wty . (72) 

— 7T 

In the limit of small "temperature" (1//3 — )■ 0) and for moderate values of m = |/c — /|, we 
use the known asymptotics of the Bessel function and obtain: 

e /3 / A m 2 — I \ 

F ^= ! vb( 1 -V + -)" (73) 

i.e., the power decay is governed by the critical exponent £ = —1/2. 

Since the summations can be replaced by the integrations at large enough M, we obtain 
from ([39]) and gU): 

pN T-rf fdd>i\ -j8E(i-«*i) 



(*,(v)i5 B ^«ff s i*,(u)> = Wn(/lr) e 



i=l 

— 7T 



(74) 

P(v- 2 , e^)P(e"^, u 2 ) H | e^ fe - | 2 , 



Kfc<KiV 



where the continuous integration variables fa £ [0, 27r] are due to the change of the discrete 
variables ki £ .M as follows: 0^ i— >■ </>j, Vi. Further, the integral f J74|) can be approximated 
at /3 tending to infinity as follows: 



(tf iv(v) | iln e"^ i7 n I V N {u)) ~ P(v- 2 , 1)P(1, U 



2\ 



e 



/3iV 



(2vr)^X! 



OO OO OO jv (75) 

X / / ' ' • / e i=1 I J | 0A; - 0i | 2 #1#2 • • • C?0JV • 



-oo— oo — oo 



Kk<KN 
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The bold-faced argument 1 in the sum, say, V(l, u 2 ) implies that S\(e ±z ^) in (j48l) is replaced 
by S\(l, 1, . . . , 1) provided the iV-tuple of the exponentials e ±%< ^ is substituted by iV-tuple 
of unities. The corresponding value of the Schur function can be obtained |56] : 

n (a, - 1 - x k + k) 

5a(M '---' 1) = (^-T)!(iV-2)!..-l!0!- (76) 

Right-hand side of ( 176]) coincides with the dimensionality = dim7r,\ of the unitary irre- 
ducible representations of the unitary group U(N), which corresponds to the signature A, 
i.e., d\ = S\(l, 1, . . . , 1) [5B]. The Schur functions provide a base of the ring of symmetric 
polynomials of N variables [56J. 

The integral in (P75i) is the Mehta integral [59] of the Gaussian Unitary Ensemble of 
random matrices. Its value is known, and the estimate for right-hand side of (T75]) acquires 
the form: 

<<Mv) | n n e-P E ** n n | * N (u)) ~ P(v- 2 , 1) V(l, u 2 ) ^ J] . (77) 

Let us consider the case when u 2 and v 2 in (T77j) are the solutions of the Bethe equations 
(flQj) . Approximate estimates u 2 = e %e ~ 1 and v -2 = e~ l ° ~ 1 are valid at 1 <C iV <C M. 
Thus, the estimate (177)) takes the limiting form: 

(v N (i, ...,i)\n n e-pxx* n n | v N (i, . . . , i)) ~ 

v h h V^LttJM. (78) 

\C{(M—N—n) N } /3^/ 2 11 ( 27 T)V 2 " 



The present estimate is proportional to the square of (jTOj) at q — 1, i.e., to the square of 
the number of boxed plane partitions of the size N x N x (M — n — 1). As a result, the 
ratio of the average T(0, n, (3) (1251) to the partiton function of the XX model Z [35] can be 
estimated as follows: 

T(0,n,p) f^r M-n + j - fc\2 e - l3£ o x 



const x f J| 



where £q * s ^ ne g roun d state energy of the XX model. In the same limit, the generating 
function (112)) is specified as follows: 



v k 



(^(1, . . . , 1) | TI n e-&™ II n | ^(1, ...,!)) 



E \P K (» R ^» L )\d xL d x n, 

\ L ,\ R C{(M-N-n) N } 

where a relationship between strict and non-strict partitions is valid: ^i L,R = X L ' R + S (see, 
for instance, fl39|) ). 

In the case of strong anisotropy, we use (1711) and obtain an analogous estimate which 
demonstrates the proportionality to the square of the number of strict plane partitions in a 
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box of the size N x N x (M - n - N): 



T(0,n,(3) 



Z 



~ const x 



( E 



\C{(M-2N-n+l) N } 




N 



M - N +l-n+j 




= const x 



(n 



j + k-i 



k,j=i 



where Z is the partition function of the model in the limit of strong anisotropy J5j[8] , and 
£q a is the corresponding ground state energy. The relation (I58p is specified as follows: 



where n L,R = * L ' R + 2 <* (see Ol ). 

6 Discussion 

The XX Z Heisenberg chain has been considered for two specific limits of the anisotropy 
parameter: A — > and A — > — oo. The corresponding state-vectors have been expressed 
by means of the symmetric Schur functions. Certain expectation values and thermal cor- 
relation functions of the ferromagnetic string operators have been calculated over the base 
of iV-particle Bethe states. The expectation values obtained are of the type of the empti- 
ness formation probability. The thermal correlator of the ferromagnetic string operator is 
expressed through the generating function of the lattice paths of random walks of vicious 
walkers. The thermal correlator in question turns out to be a generating function of certain 
polynomials build up from the Schur polynomials. A relationship between the expectation 
values obtained and the generating functions of boxed plane partitions is discussed. Asymp- 
totic estimate of the thermal expectation value of the ferromagnetic string is obtained in 
the limit of zero temperature for A = 0. These estimates are expressed in terms of the 
dimensionality of the irreducible representations of the group U(N). 



vf /2 (**(!,..., i) | n n e-^n n | **(!,..., 1)) 




p K (ji R ^ii L )\ dfrd 
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